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Abstract 

We study the complexity of approximating the smallest eigenvalue 
of a univariate Sturm-Liouville problem on a quantum computer. This 
general problem includes the special case of solving a one-dimensional 
Schròdinger equation with a given potential for the ground state en- 
ergy. 

The Sturm-Liouville problem depends on a function q, which, in 
the case of the Schròdinger equation, can be identified with the po- 
tential function V. Recently Papageorgiou and Wozniakowski proved 
that quantum computers achieve an exponential reduction in the num- 
ber of queries over the number needed in the classical worst-case and 
randomized settings for smooth functions q. Their method uses the 
(discretized) unitary propagator and arbitrary powers of it as a query 
("power queries"). They showed that the Sturm-Liouville equation can 
be solved with C(log(l/e)) power queries, while the number of queries 
in the worst-case and randomized settings on a classical computer is 
polynomial in 1/e. This proves that a quantum computer with power 
queries achieves an exponential reduction in the number of queries 
compared to a classical computer. 

In this paper we show that the number of queries in Papageorgiou's 
and Wozniakowski's algorithm is asymptotically optimal. In particu- 
lar we prové a matching lower bound of í)(log(l/e)) power queries, 
therefore showing that 0(log(l/e)) power queries are sumcient and 
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necessary. Our proof is based on a frequency analysis technique, which 
examines the probability distribution of the final state of a quantum 
algorithm and the dependence of its Fourier transform on the input. 

1 Introduction 

This paper deals with the solution of the Sturm-Liouville problem on a quan- 
tum computer. Quantum computers have shown great promise in solving 
problems as diverse as the discrete problems of searching and factoring pUÍTïï] 
and the continuous problems including integration, path integration, and ap- 
proximation |13ll5tfT6tl6ll7]. The main motivation for quantum computing is 
its potential to solve these important problems efficiently. Shor's algorithm 
achieves an exponential speedup over any known classical algorithm for fac- 
toring, but until the classical complexity of factoring is proven, the exponen- 
tial speedup remains a conjecture. The quantum algorithms for integration 
provide provable exponential speedups over classical worst-case algorithms, 
but only polynomial speedups over classical randomized algorithms. 

Recently Papageorgiou and Wozniakowski introduced a quantum algo- 
rithm for the Sturm-Liouville problem [H] which uses the quantum phase 
estimation algorithm. They showed that quantum algorithms with power 
queries 1 achieve a provable exponential reduction in the number of power 
queries over the number of queries needed in the classical worst-case or ran- 
domized setting. Naturally query complexity results neglect the cost of ac- 
tually implementing the queries. At the end of this paper we will discuss 
this problem for power queries, but it is currently not clear under which 
conditions power queries are sufficiently inexpensive to implement for the 
Sturm-Liouville problem. 

In this paper we will prové lower bounds on the number of power queries 
for quantum algorithms that solve the Sturm-Liouville problem. This can 
be used to show the optimality of the algorithm proposed in [H]. To prové 
lower bounds for algorithms with power queries the previously known quan- 
tum lower bound techniques, such as the "polynomial method" of Beals et. al 
[H [TT] do not suffice. Our lower bound method builds on the "trigonometric 
polynomial method" [2], which is an extension of the above-mentioned poly- 
nomial method and was modified to be used with power queries in [B] to prové 

J We will define power queries rigorously in Definitionn] Informally they are just an 
arbitrary (integer) power of a specific unitary matrix. 
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lower bounds for the phase estimation algorithm. Our method uses frequency 
analysis instead of a maximum degree argument, which is not applicable in 
the case of arbitrary powers. 



2 The Sturm-Liouville eigenvalue problem 

Papageorgiou and Wozniakowski study in [H] a simplified version of the 
univariate Sturm-Liouville problem. Consider the eigenvalue problem for the 
differential equation 

—u"{x) + q(x)u(x) = Xu(x) 

u(0)=«(l)=0 W 

for a given nonnegative function q belonging to the class Q defined as 

Q = íq : [0, 11 [0, 1] : q E C 2 ([0, 1]) and max max \q {i) (x)\ < l). (2) 
l 1=0,1,2 xe[o,i] J 

We are looking for the smallest eigenvalue A such that there exists a non- 
zero function u\ that satisfies (JTJ) . What is the minimal number of queries of 
q that permits the determination of the smallest eigenvalue À in this equation 
with error e and probability 3/4 on a classical or quantum computer? 

The one-dimensional time-independent Schrodinger equation 

h 2 d 2 

; ^(x) + V(x)*(x) = E*(x) (3) 



2m dx 2 



of a particle in a box, see [TH], is an instance of (JT}. We are given a potential 
V and are looking for the eigenfunctions \I/ of this equation and their corre- 
sponding energies E. In particular, we are interested in the ground-state and 
its energy, i.e., for a given potential V, we want to determine the eigenfunc- 
tion \l>o and its energy Eo, such that all other eigenfunctions \l/ n have higher 
energies E n > E Q . Since quantum systems obey equation (jHJ, it seems plau- 
sible that quantum computers could potentially solve the eigenvalue problem 
faster than a classical computer. 

In the next section we define a quantum algorithm with power queries. 
We especially have to tackle the question concerning the form of the input 
(i.e., the function q in the Sturm-Liouville problem) enters the quantum 
algorithm. 
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3 Quantum algorithms for the Sturm-Liouville 
problem 

Let us denote the differential operator associated with the Sturm-Liouville 
problem for a certain q G Q as ~h q : C 2 ([0, 1]) — ► C°([0, 1]), defined by 



h g u(x) 



dx 2 



u(x) + q(x)u{x). 



We discretize h q by approximating the second derivative at the points ^r-j-, 

M q : 



-tti • • "~ tt and obtain an n x n matrix M Q : 

n+l ' ' n+l y 



(n+l) 



2 -1 
-1 2 -1 



-1 2 -1 
-1 2 



+ 



' n— 1 N 
vn+1' 



- 



(4) 



The eigenvalues of L 9 and M q are closely related. Let us denote the smallest 
eigenvalue of h q by À(ç) and let us write \\{M q ) for the smallest eigenvalue 
of M g . Then (see e.g. |Sj) 



\{q)-\ 1 {M q ) = 0{n~ 2 ). 



(5) 



The input q G Q enters the quantum computer in the form of a unitary 
black-box transformation called a quantum query. For the Sturm-Liouville 
problem we define this query to be the unitary operator exp(|M g ). One 
can show that the smallest eigenvalue X(q) of the Sturm-Liouville equation 
satisfies tt 2 < X(q) < n 2 + 1. To avoid ambiguity we use proper scaling, i.e., 
instead of exp(iM q ) we use exp(|M g ), which defines a unique phase (p G [0, 1) 
by 2iri<p = |A(í). 

We now define an associated quantum power query for exp(|M ç ). 

Definition 1. Let L ç be the differential operator for a Sturm-Liouville prob- 
lem and M q its discretization at n points as in for q G Q. We define 
the power query Wf (exp(|M 9 )), where l G {1, 2, ... c} and p G N, acting on 

c 2C ®c r > " 



as 



W?{eMW q ))\x 1 )...\x c ) 



\xi) . . . \x c ) exp(|pM g ) |^) for x t = 1 
\xi) . . . \x c ) \ip) otherwise 
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for all xi,...,x c G {0,1} and arhitrary normalized vectors G C n and 
extend this defi,nition to all quantum states by linearity. 

Suppose that the \i/j s ), s — 1, . . . ,n, are the eigenvectors of M q and that 
M q \if) 8 ) = X s \4> s ). Then for \ip) = Y^=i a s \i*s) and \x) = \x%) . . . \x c ) with 
Xi = 1 

n 

Wf (exp(|M 9 )) \x) |V) = \x) exp(fpM g ) = k) ^ pAs |V*> • 

3=1 

Quantum algorithms are products of unitary transformations. Every 
quantum algorithm that approximates X(q) can be divided into stages that 
use powers of exp(|M ? ) and therefore depend on q, and stages that are in- 
dependent of q. Let us define a quantum algorithm with power queries. 

Definition 2. For a Sturm-Liouville problem given by the input q G Q with 
the solution X(q), we define a quantum algorithm 

A = (\^ 0) );U ,...,U T ;h, Pl ,...,l T ,p T ;X) 

with T power queries that solves this problem as follows. Let Uq, U\, . . ., Ut 
be arhitrary but fixed unitary transformations and \ip^} a fixed initial state. 
Let W[. J (exp(|M ? )) be a power query as in Definition^ A measurement of 
the state 

|^ T )(exp(fM q ))) = tf T Wf/(exp(fM g )) . . . ÍW^exp^ M q ))U \^ 0) ) 

in the Standard basis yields a state \k) with probability P k(q)- For each k 
compute an approximation X(k) G M to the eigenvalue of interest X(q) on a 
classical computer. For every q G Q the probability that an e- approximation 
X(k) of X(q) is computed is given by 

fc:|A(ç)-A(fc)|<e 

For any algorithm A with T power queries we define 

e(A,T) = inf {e : e chosen such that is larger than | for all q G Q} 
as the worst-case quantum error of A. 

We measure in the Standard basis for convenience only; a measurement in 
any other basis is easily achieved by modifying the operator Ut accordingly. 

A model like this was introduced in [T] for discrete inputs q. It was ex- 
tended to continuous functions by Heinrich in [H]. Our model is an extension 
of this model to incorporate power queries. 
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4 Upper bounds 

To estimate X(q) on a quantum computer with power queries Papageorgiou 
and Wozniakowski used the quantum phase estimation algorithm, see e.g. 
[T2] . This algorithm takes a unitary transformation Q with an eigenvector 
|£) as input, i.e., Q |£) = e 2 " v Here <p E [0, 1) is called the "phase" of the 
eigenvalue corresponding to and the phase estimation algorithm gives us 
an approximation (p to (p. This algorithm has the final state 

\^ T \Q)) = {F£ ® I)Wf-\Q)Wf- 2 (Q) ■ ■ ■ W?{Q){H® T ® 7) |0) |£> , 
and is depicted in Figure [TJ Suppose Q is a r qubit transformation. A 

h) 

k 2 ) 

k T -i) 
k T ) 



Figure 1: The quantum phase estimation algorithm. T~t is the inverse 
quantum Fourier transform on T qubits. 

measurement of \^ T '(Q)) returns a state 

\k) = \h) . . . \k T ) \k T +i) ■ ■ ■ \k T + r ) ■ 

The algorithm then uses k to compute an approximation <p(k) = k\l~ x + 
/c 2 2 — 2 + . . . + kT2~ T to <p classically. 

One can show, see e.g. [TJ], that with probability greater than | the al- 
gorithm approximates <p up to precision e with 0(log((l/e))) power queries. 
Papageorgiou and Wozniakowski use this algorithm to approximate the small- 
est eigenvalue X(q) of the Sturm-Liouville operator L g and use the operator 
Q = exp(|M g ) as a query. Since the phases of exp(|M 9 ) and exp(|L g ) are 
related through equation (JHJ), we have to discretize at n — 0{e~ 1 / 2 ) points. 



0) -\_H_ 
0) -\H 



0) - H 



0) - H 
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The quantum phase estimation algorithm requires the knowledge of the 
eigenvector for which the phase is estimated. For the Sturm-Liouville prob- 
lem we need the eigenvector \zi(M q )) of M q corresponding to the small- 
est eigenvalue Ài(M g ). We can compute \zi(M q )) through the method of 
Jaksch and Papageorgiou [8j, which computes a superposition of eigenvec- 
tors \zj(M q )) of M q , with a sufficiently large \z\(M q )) component. For details 
see làíl- ' 



5 Lower Bounds 

Our goal is to prové that the algorithm described in the previous section 
is optimal with respect to the number of power quedes. We have to prové 
that every quantum algorithm A with T power queries that returns a correct 
answer with precision e(A, T) < e has to use T = í7(log(l/e)) power queries. 

We will show that even for a much simplified version of the problem this 
lower bound still holds. Consider as input only constant functions q(x) = q G 
[0, 1]. Obviously q G Q. It is easy to see that in this case the eigenfunctions 
which fulfill the boundary condition in JTJ) are 

u s (x) = sin(s7rx) (7) 

for s G N and that they have eigenvalues À s = s 2 tt 2 + q, which means that 
the smallest eigenvalue X(q) is A(ç) = n 2 + q. 

Similarly for the discretization M q of L g with constant q G [0, 1] the 
eigenvectors are 

n 

!«•>= v^5> n ^ |x) (8) 

x=l 

with eigenvalues 4(n + l) 2 sin 2 ( 2 ,^_ 1 \ ) + q- 

We want to investigate how different power queries lead to different out- 
puts and turn to the techniques in j3j. 

Theorem 3. Any quantum algorithm with power queries (exp(^M q )) for 
q(x) — q G [0, 1), see Definition\^ that uses c G N control qubits, can be 
written as 

\^ T \exp(iM q ))) = U T Wf T T (exj>^M q )) . . . UiW^(exp(íM q ))U \^) 

= E s í T, te) I*) ■ ( ' 

k=0 
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where U\, . . . ,Ut are unitary operators and the Sf. (q) are trigonometric 
polynomials of the following form: 

SP(q) = £ V ( k T le^, (10) 

meMT 

with ÀAt defined as M.q = {0} and 

Mt+i = {m : m e Mt} U {m + p T+ i : m G Mt) , (H) 

(T) 

and the coefficients ïj^GC do not depend on q and are normalized: 

E E 1C1 2 = i- (12) 

k m&Mr 

Proof. The proof is by induction on the number of quedes T. We will 
write the state of the algorithm after T steps \ip( T ' (exTp(^M q ))} in the basis 
(\k) \ift s ))k,si k — 0, 1, . . . , 2 C — 1, s — 1,2, ... ,n, which is split into a control 
part \k) and an eigenvector part \ip s ). We will not address the ancilla qubits 
in our proof, but they can easily be treated (after possibly reordering the 
qubits) as control bits that are never used. 
For T = power queries we can write 

|^°)(exp(fM 9 ))) = U |V (0) ) = EÍ,o IV>.) , 

fe,s 

which contains only powers e^ mq from m E Aio = {0} and obviously 

l^l«,ml = E/ l^fe,a,ol = I' 
fc,s mdMo k,s 

Let us now assume |^ T )(exp(|M 9 ))) can be written as 

l^ T Wf M q ))) = e E C!m e|m9 i^) , 

k,s mÇiMj· 

with coefficients ?7^ m fulfilling condition (JT2J). If we apply W / z ï ^ 1 (exp(|Mg)) 
to |^A T )(exp(~M g ))) we get (ki T+1 is the control bit, i.e., the Ix+i-th bit in 
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the binary representation of k): 

<- 1 (exp(fM (? ))|^)(exp(|M g ))>= £ £ V^e^ \k) \^ s ) 

k,s radMT 



k T+1 =o 



+ E E I*> exp(^T + iM 3 ) i^} . (13) 



k,s vh^Mt 

V+i =1 



We define ( s := e 2 ^"' 4 " 1 -' &m (2(71+1)) anc j p r0 ceed to analyze the second term 
in (fT3|) . where the control bit fc ÍT+1 = 1 and get the following 

E ^l m e^ mq \k)exp(lp T+1 M g )\^) 

m£MT 

= E C! m ^ m9 e^ +i ( 4(n+1)2sm2 fe) +9 )|fc)i^) 
= E CUr +lef (m+PT+l)9 • 

If we define vies'm f° r a ^ m e -A^t+i as 

{V ( k T l m - PT+1 (s T+1 for ^t+i = 1 and m - p r+i e M T 
VkJ, m for K+i =0 and me M T 

otherwise 

we can write 

<;-(exp(fM (? ))|^ T )(exp(fM 9 ))> = ^ £ l*> IV>.> ■ 

k,s mSA^T+l 

We check our normalization condition (|T2| for rf^s^m > 

k,s mGMT+i 

= E E iCJ 2 + E E l-C^„cr«l 2 

fci-r , 1 =0 h , , =1 



'T+l 



E E l^.s.ml ~~ l- 
k,s m&MT 
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The next step in the algorithm is to apply the unitary transformation 
Ut+i- For k,l = 0, . . . , 2 C — 1 and s, t = 1, . . . , n define the coefficients 

Ui,t,k,s = (l\ (-0*1 Ut+i \k) \ip s ) and let 



'll,t,m 



E-í 



(T+l) 



fe,s 



This allows us to write 

^ +1 < T ^ 1 (exp(fM g ))|^( T )(exp(fM (? ))> 

=E E €í? e|m? ^ + i i*> 



E E EC^W> 9 ioi^> 



(T+l) 



E E 

l,t meMT+i 



ll iq 



It remains to check that 



E E 

Z,í mgA^T+i 

E E 



L(t+i) 



k,s 



UL t ,k,. 



ec; 1 ' 



Ult,k',s' 



k's' 



E E 

k,s,k',s' meMT+i 

E E 1^ 



<T+1) 
lk,s,m 



. z,t 



~{T+1) 
^lk' ,s' ,m 



{T+l) 
s,m 



where we used that C/t+i is unitary. This completes the proof. 



□ 



We can use Theorem to get explicit formulas for the probability of 
measuring a certain state. 

Lemma 4. Let A be a T power query quantum algorithm for the Sturm- 
Liouviüe problem with powers pi,...,px and c G N control bits as defined in 
Definition\^ Let B be a partition of the set of all basis vectors, i.e. 

(J B = {\k) : k = 0,1,..., n2 c - 1} and BnC = for B,C G B, B^C. 
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// the input q G Q is a constant function q(x) — q G [0, 1), the probability of 
measuring a state \k) from B G B is a trigonometric polynomial 



iec T 



(14) 



with coefficients f3 B ( G C that are bounded by 

E K/i si 

BeB 

for all possible partitions B, and the set Ct is given by Cq = {0} and 

C T +i= (J {1,1+Pt+i,1-Pt+i}- 
leC T 



(15) 



Proof. Consider quantum queries exp(|M g ) for constant functions q(x) = 
q G [0, 1) in the Sturm-Liouville problem. From equations ((HI), (fTÜ|) we know 
that the final state of every T power query algorithm can be written as 

k meMr 

Let B be a partition of the set of all basis states \k). Thus the probability to 
measure a state from the set B G B is 

2 



Mg) = E 

keB 



E C 



,e2' 



feeB 



E (C)* 6 " 1 "" 1 



E -í 

.m2eA^T 



( T ) e 2 m2 9 



iri2 



E E l'/À'»'nj 'li,,,/ 2 

k€B m 1 ,m2€MT 



( T ) „ïï(m2-mi)9 



with coefficients / defined as 



3 {T) 



E E 

fcGB mi,m2£MT 
1712— m\=l 



ST) 



(T) 



(16) 
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and the set Ct is given by 

Ct = {mi — m 2 : mi, m 2 G Mt} ■ 

(T) 

For any partition B we can now bound the j3 B ( as follows 



(17) 



EK?I=E E E 



BeB 



BeB k&B m 1 ,rn 2 £MT 
rri2—mi=l 



Vic,™,! j Vk 



:,m 2 



sE E 



fe mi,m2eAlT 
m2-rai=l 



where Y^ fc is the sum over all possible states \k). From (fT2|) we now derive 
by the Cauchy-Schwarz inequality 



E E = E E 

k m:m,m+l&MT 



fe mi,m2GAlT 
mi— mi=l 



(T) (T) 

,771+2 



1/2 



< 



E E h£ 

fe VmeA'ÍT 



(T) 
m 



m+Ze.MT 



1/2 



£ E E K 

fe meAlr 



(T) 
rei 



< 1. 



It remains to show that the two definitions of Ct in equations (fT5|) and (|T7l) 
are identical. The proof is by induction. T = is trivially true. We use the 
definition (fTT]) of Mt to see that 

C T+ i = {mi - m 2 : m ll m 2 G M T +i} 

= {nii - m 2 , mi + p T+ i - m 2 , m x - m 2 - Pt+i, 
mi + p T +i -m 2 - p T+ i : mi, m 2 G Mt} 
= {l,l + Pt+i, l - Pt+i ■ l G C T } , 



which completes the proof. 



□ 



Note that \Cr\ < 3 T . This bound is sharp, since for the choice of pi = 3* 1 
we have C = {0}, A = {-1,0,1}, C 2 = {-4, -3, -2, . . . , 3, 4} and in 
general 



{-3 T - 3 T_1 - . . . - 1, . . . , 3 r + 3 T_1 + . . . + 1} . 
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5.1 Fourier Analysis of Power Query Algorithms 

With Theorem El and Lemma|3]we have the tools needed to provide a lower 
bound for the Sturm-Liouville problem. We are now able to apply our fre- 
quency analysis technique to this problem. 

Theorem 5. Any quantum algorithm A with T power queries which esti- 
mates the smallest eigenvalue X(q) in the Sturm-Liouville eigenvalue problem 
for all inputs q(x) = q G [0, 1) with precision e(A,T) < e and probability 
greater than 3/4 has to use T = í7(log(l/e)) power queries. 

Notice that a lower bound on the "easy" subset of constant functions 
q(x) = q implies that the same lower bound holds for any set of inputs that 
includes the constant functions, hence it also holds for the class Q. We also 
would like to remark that the lower bound T = Q(\og(l / e)) does not depend 
on the number of discretization points n. 

Proof. After T power queries we measure the final state and receive a state 
\k) with probability Pk(q)- From the integer k we classically compute a 
solution X(k). A successful algorithm has to return an e-approximation for 
every q G [0, 1) with probability 

S Pk ^ - 4' 

k:\\(q)-\(k)\<e 

see Definition [2J Define 

A q , £ :={k: \X(q)-X(k)\ < e} 

as the set of states that are mapped to e-correct answers for input q. Choose 
iV G N such that is slightly bigger than 2e, i.e., jÀ-y < 2e < ^ and define 
the points x r := (r + 1/2) /N for r — 0, 1, . . . , N — 1. For the inputs q = x r 
we can visualize the quantum algorithm A as in Figure EJ Notice that the 
sets A Xrt£ are mutually disjoint for r = 0, . . . , N — 1, because x r and x r+ \ are 
chosen such that 



|A(ar r ) - A(x P+1 )| = 16sin 2 (f ) + r -jf- - 16sin 2 (f) 



N 



Therefore there can be no state \k) that is mapped to an output A(fc), which 
is an e-approximation to X(x r ) and X(x r+ i) at the same time. Let 



13 



Input q 




Output: approximation A to >^(q) 



Figure 2: A quantum algorithm for the Sturm-Liouville problem with inputs 
q — x r , r — 0, . . . , N — 1, will result in a probability distribution Pk(q) on the 
states \k) that are measured. Each state \k) is mapped to an answer \{k). We 
write A Xrjt for the set of all states \k) that are mapped to e-approximations 
of \(x r ). 

be the probability of measuring an e-approximation to X(x r ). Since the sets 
A Xry£ partition the set of all outputs, Lemma |U allows us to write 

iec T 

We apply the iV-point inverse discrete Fourier Transform to p r>£ (ç), which 
we evaluate at the points x n , and get the following value at k = 0, 1, . . . , N—l: 
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N-l 



DFT N [p rt£ }(k) =J2 



■2-!rikn/N 



n=0 
N-l 



££<C! 

n=0 1&C T 



(T) e ~l{n+l/2)/N o -2nikn/N 



N-l 



(19) 



2ni( 1 L· -k)n/N 



n=0 



2-irilJ k) ^ , 



iV 



^ (mod JV) 
,^ = À;(mod iV) 



where = k (mod iV) indicates that there exists an integer z such that 
= k + zN. For every l define l/4 n rm G [0, N) as 

l /MN) := min {//(4tt) - zN : z = 0, 1, 2, and //(4tt) - ziV > 0} . 

Then exp(27rz^/iV) = exp(2Tcil/ Í7r ^ N )/N). To take the absolute value of 
equation (fT9f . we use |e í2ÉI — 1=2 |sin(0)| and get 



\DFT N \p r)£ ){k)\ < 



£ UT) 

ie£ T 



3 K ' 



\sm(ir(l/4ir(N)-k))\ 
| sïn(7r(Z /47r(JV) — fc)/A^ ) | 

N 



l/4ir(N) 7^ k 
l/4n(N) — k 



(20) 



We can bound the Fourier transform (fTÏÏ|) by separating the correct an- 
swers, i.e., the e-approximations to x r , from the rest: if the input q = x r 
then the algorithm has to return an answer À that is e-close to the correct 
answer X(x r ) with probability greater or equal 3/4. This probability is given 
by p r ,e{q), i-6., we demand that p re {x r ) — 3/4. Then: 



N-l 
n=0 



-2-Kikn/N 



> 



Pr,e (*^r ; 



JV-1 
n=0 



>J - J2PrA X n), 
n=0 



N-l 



(21) 
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Consider the second term in f2*Tl) . ^2n=o Pr,e( x n) ■ Recali that p r ,e(l) is the 

probability that the algorithm measures a state \k) that is mapped to an 
answer X{k) that is an e-approximation to X(x r ), i.e., \k) e A Xr>e , see (fT8j) . 
This probability p r ,e(q) depends on the actual input q. For input q = x n ^ x r 
a state \k) G Ac^s will noi yield an e-correct answer: we chose the x n , 
n = 0, . . . , N — 1, such that \X(x n ) — A(x r )| > 2e for n/r, and thus there 
cannot be an e close answer for both x r and x n . The sum Yun=o Pr,e{%n) now 

telis us how often the algorithm chooses a state from A TjE . 

If we knew that none of the wrong answers is preferred by our algorithm, 
say e.g. Y^n=o Pr,s{x n ) < \, equation ([2Ï]) would read 

n^ír 

, N-l N-l 

E, , 2-KÍkn O V' "\ , < 1 / \ 

PrA X n)e N > J ~ 2_^PrA X n) > (22) 
n=0 n=0 

We will show that this property has to be true for some r = 0, . . . , N — 1, 
indexing the set of states A Xr>£ that represents numbers e-close to x r . Let 
i? < be the set of all r for which J2n=o Pr,e{ x n) < \ holds and i?- the set for 

which it does not. We estimate the number of elements of i? < by splitting 

N-l N-l N-l N-l N-l N-l 

N = Y 1 - Yl Y pr ^ Xn "> = 5Z pr ' £ ( Xr ) + YY pr ^ x ^ 

n=0 n=0 r=0 r=0 r=0 n=0 

into the following parts: 

N-l N-l N-l 

N > ^Pr, £ (x r ) + Y Y Pr ' £ ( Xn "> + zJ Y Pr > £ ( Xn "> 
r=0 rgi?.< n=0 rG-R- n=0 

>N- +| J R < |·0 +1^-1- 

4 1 1 2 

and therefore we can conclude that < \N and thus \R<\ > \N . Now 
> implies that we can actually choose an element r e i? < . Fix such 
an r and we can combine equations (|2Ü|) and (|2*2*|) to 



sin( 7r ( í /47r(JV)-fc)) , , 



V4<E K3 i l·l^WF^ 'W^* . (23) 
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We will now fix the parameter k = 0, 1, . . . , iV — 1 in inequality (Í2~ÏÏ1) in such a 
way that the terms in the sum on the right-hand-side of the inequality are as 
small as possible. This will imply that the sum must be over a large number 
of elements, i.e., that \Cr\ is large. Since \Cr\ < 3 T this will help us to 
ultimately prové that T = í2(logiV) if we could show that there is an a > 
such that \Ct\° = íí(iV). More specifically we will show that |£t| 2 > 
which proves T = f2(logiV). 

l 1 2 1 l 1 2 1 

We prové \Ct\ > j^N by contradiction. Assume \Ct\ < JqN. This 
assumption allows us to find a k such that the right-hand-side of inequal- 
ity (|23|) is smaller than the left-hand-side, which will lead to our desired 
contradiction. 

If we project Ct into the interval [0, N) through l \— > I/^n) we will get 
a set {1/a-k(n) ■ l G £t}- Order this set as < t\ < í 2 < • • • < t\c T \ < ^· 
This defines "gaps" between these numbers, i.e., intervals G = (tj, íj+i) for 
j — 1, . . . , \Ct\ if we define t\c T \+i = íi + N (we "wrap around"). Define the 
width w(G) of such a gap G as the distance between its endpoints. Thus 
w({tj, t j+1 )) =t j+1 -tj. 

Let G m be the gap with the maximal width w{G m ) in the distribution. 
Its width must be w(G m ) > N/ \£t\, since 

N = J2 W ( G ) < ^ max ^( G ) = \C T \maxw(G). 

G G 

Additionally w(G m ) > 10, since we assumed |£r| 2 < jqN and therefore 
t^t > 10 \Ct\ > 10. Thus there are at least ten integers k G {0, 1, . . . , N — 1} 
that fall into this largest gap G m , i.e k G G m . One of these k has maximum 
distance to both boundaries tj and t,+i of G m : it is the k that is closest to 
the middle m = tj+1 ^ tj of G m = (tj, tj+i). This integer k fulfills \k — m\ < \ 
and 

min {k — tj, tj + i — k} = min {m — tj + k — m,m — k + tj + \ — m} 

> w(G m ) 1 > N 1 
2 2 ~ 2\Ct\ 2' 

Fix this G (íj, íj+i). Now |sin(x)| > 2/7r|x| for —tt/2 < x < n/2 and 
therefore 



mm 

iec T 



. n(li ilT (N) - k) 
sm — 



N 



2 i i 2 

> min — \l/4n(N) -k\ = — min {k - tj, t j+1 - k} 

1 1 

-LC^ ~ ÏV' 
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Then we can use this to estimate (|2~ÏÏ|) : 



1/4 < E |/C 



|sin(7r(// 47r (jv) - fc))| 
|sin(7r(Z /47r(A r) - fc)/iV)| 



£ E C 



1/|£t| - 1/iV 



We sum the last inequality over all r G R < for which it is vàlid, and get: 

1 



4 - i 



reR< 



/\Ct\- l/N 



E E \fi 

reR< 1&C T 



(T) 
r.s.l 



(24) 



Since the number of elements in R K is bounded by \R < \ > \N , the left-hand- 
side of ((2H) is bounded by |iV < \R < \ \. The right-hand-side of inequality 
JZH> can De bounded through LemmalU 



E E |a 



(T) 
r,e,í 



< \Cn 



1 2 i 

If we put both sides together again and recali that we assumed \Ct\ < jqN 
we get 



8 ~ 











2 


l/l£ï 




■l/N 1- 




\/N 



< 



J-jV 

10 



< _ io 



10\Cj 



i 9 

10 



In, 



which is a contradiction. 

Therefore \Ct\ 2 > jqN must hold. This, together with \Ct\ < 3 T , leads 
us to iV < 10 • 9 T . Take the logarithm and we get T = f2(logiV). We chose 
N such that < 2e < which finally proves that the number of power 
queries T for any algorithm A with error e(A, T) < e has to be of the order 
r = íí(log(l/e)). □ 



6 Discussion 

In this paper we have proven lower bounds for the number of quantum power 
queries for the Sturm-Liouville problem and settled an open problem in [14J. 

How does this number of T = 0(log(l/e)) power queries relate to the 
cost of quantum algorithms? Here we understand "cost" as an abstraction 
on the number of elementary quantum gates or the duration for which a 
Hamiltonian has to be applied to a quantum system. Suppose the function 
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q is from a class Q' Ç Q where each power query Wf (exp(~M ç )) can be 
implemented with cost(Wf(exp(|M g ))) = cost(Q',p). 
If we implement VPjf(exp(|M g )) naively as 

^(exp(|M,)) = (^ 1 (exp(fM g ))) P , 

then cost(Q',p) = p ■ cost(Q') and the cost of the Sturm-Liouville algorithm 
with T = 0(log(l/e)) power queries grows as 

T-1 T-1 

cost(Q', 2 j ) = 2i • cost(Q') = (2 T - 1) • cost(Q') = 0(l/e) • cost(Q'). 

This is polynomial in 1/e just like the Sturm-Liouville algorithm with bit 
queries discussed in [H]. To take advantage of the proposed power query 
algorithm it is therefore necessary to realize power queries Wf(exp(iM q )) on 
a quantum computer in such a way that cost(Q',p) = o(p) ■ cost(Q') 

The implementation of power queries with cost that is not linear in the 
power p of the query is still not settled and requires more work. It would 
be of interest to identify subclasses Q' Ç Q for which we are able to prové 
cost(Q',p) = o(p) • cost(Q'). 

Another open question is whether it is possible to extend the methods 
we used for upper and lower bounds for the Sturm-Liouville problem in one 
dimension to similar problems in higher dimensions. Most important for this 
problem is probably the extension of the results in |H] on approximations of 
the eigenvector with the smallest eigenvalue to higher dimensions. 
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